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Abstract 

We investigate lepton flavor violating decays in a SUSY SO(IO) model with symmetric textures recently 
constructed by us. Unlike the models with lop-sided textures which give rise to a large decay rate for /i — > 67, 
the decay rate we get is much suppressed and yet it is large enough to be accessible to the next generation 
of experiments. We have also investigated the possibility of baryogenesis resulting from soft leptogenesis. 



We find that with the soft SUSY masses assuming their natural values, B' = \/BMi ^ 1.4 TeV and 
Im{A) ^ 1 TeV , the observed baryon asymmetry in the Universe can be accommodated in our model. We 
have also updated the predictions of our model for the masses, mixing angles and CP violating measures in 
both charged fermion and neutrino sectors, using the most up-to-date experimental data as input. 
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I. INTRODUCTION 



After Neutrino 2004, the allowed region for the neutrino oscillation parameters has been reduced 
significantly, and their measurements have now entered the precision phase. There have been a few 
supersymmetric (SUSY) S'0( 10) models constructed aiming to accommodate the observed neutrino 
masses and mixing angles (For a recent review on SO(IO) models, see Ref. jll.) By far, the LMA 
solution is the most difficult to obtain. Most of the models in the literature assume "lopsided" 
mass matrices. In our model based on SUSY S'O(IO) x SU{2) P] (referred to "CM" herein), we 
consider symmetric mass matrices which result from the left-right symmetric breaking of 5*0(10) 
and the breaking of family symmetry SU{2). In view of the much improved experimental data on 
neutrino oscillation parameters as well as those in the quark mixing from B Physics, we re- analyze 
our model and find that it can still accommodate all experimental data within la. We investigate 
several lepton flavor violating (LFV) processes in our model, including the decay of the muon 
into an electron and a photon, which is the most stringently constrained LFV process. We also 
investigate in this paper the possibility of baryogenesis utilizing soft leptogenesis. 

This paper is organized as follows: In Sec. ^ we briefly describe our model, and show its 
predictions for the masses, mixing angles and CP violating phases in both charged fermion and 
neutrino sectors, using the most up-to-date experimental data as input. Various decay rates for 
lepton flavor violation processes are calculated in Sec. II I II Sec. IIVI concerns soft leptogenesis in 
our model, while Sec. ^ concludes this paper. 



II. THE MODEL 

The details of our model based on S'O(IO) x SU{2)p are contained in CM Pj.The following 
is an outline of its salient features. In order to specify the superpotential uniquely, we invoke 
.^2 X -^2 X Z2 discrete symmetry. The matter fields are 

V'a- (16,2)-++ (a = 1,2), ^3 ~ (16,1)+++ 

where a = 1,2 and the subscripts refer to family indices; the superscripts +/— refer to (.^2)^ 
charges. The Higgs fields which break S'O(IO) and give rise to mass matrices upon acquiring 
VEV's are 

(10,1): T+++, r2-+-, r3--+, T^— , T+- 

(T26,l): C , 0++^, Ct^- 

2 



Higgs representations 10 and 126 give rise to Yukawa couplings to the matter fields which are 
symmetric under the interchange of family indices. SO (10) is broken through the left-right sym- 
metry breaking chain, and symmetric mass matrices arise. The SU{2) family symmetry |^ is 
broken in two steps and the mass hierarchy is produced using the Froggatt-Nielsen mechanism: 
SU{2) U{1) nothing where M is the UV-cutoff of the effective theory above which the 
family symmetry is exact, and eM and e M are the VEV's accompanying the flavon fields given 

by 



(1.2) : 

(1.3) : 



s, 



(2) ' 
(2) ' 



The various aspects of VEV's of Higgs and flavon fields are given in CM. 
The superpotential of our model is 
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W = Woirac + 



fRR 
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ac = IpSlpsTl + J^i'si'a (72</>(l) + r3'A(2)) 

1 — 1 

+ J^i^ai^b (74 + C) 5(2) + —tpa'4'bT5S^i) 
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The mass matrices then can be read from the superpotential to be 
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where Mu = (lO+>, Md = (10^), = (lO+>/(lO+>, = (lO+>/(lO+> and p 
^126 ^/(|lOj^). The right-handed neutrino mass matrix is 
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with Mr = l\2Qi\. Here the superscripts +/ — /O refer to the sign of the hypercharge. It is to 



be noted that there is a factor of —3 difference between the (22) elements of mass matrices M^, 
and Me. This is due to the CG coefficients associated with 126; as a consequence, we obtain the 
phenomenologically viable Georgi-Jarlskog relation. We then parameterize the Yukawa matrices 
as follows, after removing all the non-physical phases by rephasing various matter fields: 

/o a^ 

he'^ c 
a c 1 

^ ee-'^ o\ 

ee^€ (1,-3)/ 
y 1 

We use the following as inputs at Mz = 91.187 GeV 



\ 



d 



(7) 



d,e 



h 



(8) 



ruu 
■rric 
rrit 
rue 
rrifj, 
rrir 

\Vus\ 
\Vub\ 

\Vcb\ 



2.21 MeV{2.33toil) 
682 MeV{677tlt) 
181 Gey(18lil3) 
0.486 Mey(0.486847) 
103 AfeF(102.75) 
1.74 Gey (1.7467) 
0.225(0.221 - 0.227) 
0.00368(0.0029 - 0.0045) 
0.0392(0.039 - 0.044) 



where the values extrapolated from experimental data are given inside the parentheses. Note that 
the masses given above are defined in the modified minimal subtraction (MS) scheme and are 
evaluated at Mz- These values correspond to the following set of input parameters at the GUT 
scale, Mgut = 1-03 x lO^^ GeV: 

a = 0.00250, b = 3.26 x lO'^ 
c = 0.0346, d = 0.650 
e = 0.74 
e = 4.036 xlO-^ / = 0.0195 
h = 0.06878, ^ = -1.52 

91=92 = 93 = 0.746 (9) 

the one-loop renormalization group equations for the MSSM spectrum with three right-handed 
neutrinos are solved numerically down to the effective right-handed neutrino mass scale, Mr. 
At Mr, the seesaw mechanism is implemented. With the constraints \m,y^\ S> Irriiy^l^ I'^i^il ^^'^ 
maximal mixing in the atmospheric sector, the up- type mass texture leads us to choose the following 
effective neutrino mass matrix 

/ t \ 

1 1 + r 
i+r 1 J 

with n = 1.15, and from the seesaw formula we obtain 

= 7 (11) 

62 = (12) 

r 

-a{be''^{l + t^-^^)-c) + bcte''^ 
S3 = , (13) 

where r = (c^t + a^t°-^^(2 + 1^-^^) - 2a(-l + c + ct^-^^)). We then solve the two-loop RGE's for the 
MSSM spectrum down to the SUSY breaking scale, taken to be mt{mt) = 176.4 GeV, and then 
the SM RGE's from mt{mt) to the weak scale, Mz- We assume that tan/3 = Vu/v^, = 10, with 
vl + vl = (246/\/2 GeV)'^. At the weak scale Mz, the predictions for aj = /47r are 

ai = 0.01663, a2 = 0.03374, = 0.1242 . 

These values compare very well with the values extrapolated to Mz from the experimental data, 
(cKi, 02,03) = (0.01696,0.03371,0.1214 ± 0.0031). The predictions at the weak scale Mz for the 



TABLE I: The predictions for the charged fermion masses, the CKM matrix elements and the CP violation 
measures. 
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25 
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93.4ti3^Mey 


m.OMeV 


rrib 


3.00 ± 0.1 IGeF 


3mGeV 


\Vud\ 


0.9739- 0.9751 


0.974 


\Vcd\ 


0.221 - 0.227 


0.225 


\VcA 


0.9730- 0.9744 


0.973 


\Vtd\ 


0.0048- 0.014 


0.00801 


\Vts\ 


0.037-0.043 


0.0386 


\Vtb\ 


0.9990 - 0.9992 


0.999 


T'' 

'^CP 


(2.88 ±0.33) X 10-5 


2.87 X lO-'^ 


sin 2a 


-0.16 ±0.26 


-0.048 


sin 2/3 


0.736 ±0.049 


0.740 


7 


60° ± 14° 


64° 


P 


0.20 ±0.09 


0.173 


V 


0.33 ±0.05 


0.366 



charged fermion masses, CKM matrix elements and strengths of CP violation, are summarized 
in Table. The predictions of our model in this updated fit are in good agreement with all 
experimental data within la, including much improved measurements in B Physics that give rise 
to precise values for the CKM matrix elements and for the unitarity triangle P|. Note that we 
have taken the SUSY threshold correction to to be —18% 0|. 

The allowed region for the neutrino oscillation parameters has been reduced significantly after 
Neutrino 2004. In the atmospheric sector, the global analysis including the most recent K2K result 
yields, at 90% CL 

Aml^ = 2.310:1 X (14) 
sin2 29atn, > 0.9 (15) 
(best fit value: sin^ 26*0^^ = 1-0) . (16) 

In the solar sector, the global analysis with SNO and most recent KamLAND data yields, at 
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la (3(7) 



tan"' Op 



i.2 



-0.3/ 
-0.3'. 



-1.0\ 
-O.8.' 







0.39 



+0.05('+0.19n 
-O.O4Uo.llJ 



(17) 
(18) 



Combining with the CHOOZ result, a global analysis shows that the angle ^13 is constrained to 
be 3 



sin^ 6^13 < 0.015(0.048) 



at la {3a). Using the mass square difference in the atmospheric sector Am 



2 

atm 



(19) 



2.33 X 10-3 eV 



and the mass square difference for the LMA solution Attiq = 8.14 x 10 ^ eV"^ as input pa- 
rameters, we determine t = 0.344 and Mr = 6.97 x W^^GeV, which yield {51,62,63) = 
(0.00120, 0.000703e* (^-^^^ 0.0210e* (O-i^s)), We obtain the following predictions in the neutrino sec- 
tor: The three mass eigenvalues are give by 



{mu,,m^^,m,y^) = (0.00262,0.00939,0.0492) eV 



(20) 



The prediction for the MNS matrix is 



\Umns\ 



852 0.511 0.116 



0.427 0.560 0.710 
0.304 0.652 0.695 



(21) 



which translates into the mixing angles in the atmospheric, solar and reactor sectors. 



sin^26'atm = 



tan^ 9(7) 



sm ^13 



4\u Pir/ P 

\U P^2 



(1 



1.00 



I'^eui I 



\u, 



eu3 I 



0.36 



0.0134 . 



The prediction of our model for the strengths of CP violation in the lepton sector are 

fcp = Im{UuU^2U2iU22} = -0.00941 
(a3i,a2i) = (0.934,-1.49) . 



(22) 

(23) 
(24) 



(25) 
(26) 



Using the predictions for the neutrino masses, mixing angles and the two Majorana phases, 031 
and 021, the matrix element for the neutrinoless double P decay can be calculated and is given by 
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|<?7i>|=3.1xlO ^ eV, with the present experimental upper bound being 0.35 eV [4]. Masses 
of the heavy right-handed neutrinos are 

Ml = 1.09 X 10^ GeV (27) 
Ms = 4.53 X 10^ GeV (28) 
Ms = 6.97 X 10^2 GeV . (29) 

The prediction for the sin^ ^13 value is 0.0134, in agreement with the current bound 0.015 at la. 
Because our prediction for sin^ ^13 is very close to the present sensitivity of the experiment, the 
validity of our model can be tested in the foreseeable future [lo| . 

III. LEPTON FLAVOR VIOLATING DECAYS 

In light of the neutrino oscillation, extensive searches for lepton flavor violation processes, such 
as £i — > £,'7, i~ iJ^t^J: muon-electron conversion, are underway. In the SM, as the lepton 
number is conserved, there is no lepton flavor violation. Non-zero neutrino masses imply lepton 
number violation. If neutrino masses are induced by the seesaw mechanism, new Yukawa coupling 
involving the RH neutrinos can induce flavor violation ll|, similar to its quark counter part. In 



the non-supersymmetric case, the decay amplitudes for these processes are inversely proportional 
to the RH neutrino mass, M^, which is typically much higher than the electroweak scale. As a 
consequence, in non-supersymmetric models, these processes are highly suppressed to the level that 
are unobservable. 

Significant enhancement in the decay rate can be obtained in supersymmetric models, as the 
characteristic scale in this case is the SUSY scale, which is expected to be not too far from the 
electroweak scale. Thus the amplitudes for these decay processes scale as inverse square of the 
SUSY breaking scale, rather than 1/M^. The relevant interactions that give rise to lepton flavor 
violating decays come from the soft-SUSY breaking Lagrangian, 



+J^iHd^i.eL^ + ^B'JuRMR^ + BhHdHu 
+h.c. 



(30) 



where II, and vr are the LH slepton doublets, RH charged sleptons, and RH sneutrinos, 
respectively; {Hu) and (Hd) are the two Higgs (higgsino) doublets in MSSM. Assuming 



mSUGRA boundary conditions at the GUT scale, 



imj)ij = im~)ij = {m~)ij = mo6ij (31) 



= = ^0 (32) 



Ai^ = {Y,)^jAo, = {Y,)i,Ao (33) 
B'J = M.^^Bo, Bh = tiBo (34) 

where Yi, and Yg are the Yukawa couphngs of the neutrinos and charged leptons, and My^^ is the 
Majorana mass matrix of the RH neutrinos. As the slepton mass matrix {rn'~)ij is flavor-bhnd at 
the GUT scale, there is no flavor violation at Mqut- However, as {Tn'~)ij evolves from Mqut to 
the RH neutrino mass scale, M/j, according to the renormalization group equation, 



mliY^Y^h 



dln^i' L^'' 16^2 

+2(yJmlY,)ij + ml(YjY,)ij 

+(A'lA^)ij 



for i / j , (35) 

the off diagonal elements in the slepton mass matrix m~ can be generated at low energies due to 
the RG corrections jl^ . 



Svr 



X iylMyu)kjH^^) , (36) 



k=l,2,3 

for i ^ j- Here is the Yukawa couplings for the neutrinos in the basis where both charged lepton 
Yukawa matrix and the Majorana mass matrix for the RH neutrinos are diagonal; M/j^ are the 
masses of the heavy neutrinos. The Yukawa coupling in the new basis is related to Y^ in the 
original basis by 

y. = PrOrY.oI^. (37) 
Here Ol^ is the diagonalization matrix for 

Mf^^ = Oe^MeOl^ , (38) 

and the diagonal phase matrix Pr and the orthogonal matrix Or are defined by, 

Mt;;% = diag(Mi,M2,M3) 

= PrOrM.^^OIPr , (39) 



TABLE II: Summary of current status and future proposals of the experimental searches for lepton flavor 
violating decays. 



Decay 



current bound on the branching ratio 



reach of future experiment 



3e 



II e in 22^* 



T ■ 

r ■ 



67 



< 1.2 X 10"" (MEGA, 1999)[13j 

< 1.0 X 10-12 (SINDRUM, 1988)[1^ 

< 6.1 X 10-13 (SINDRUM II, 1998)[16] 

< 3.1 X 10-^ (BELLE, 2003) [18] 

< 3.6 X lO-'^ (BELLE, 2003) [19] 



10-" (PSI)[14] 
10-15 (J-PARC) 

2.0 X 10-1^ (MEC0)[17] 
10-i« (J-PARC) 
10-9 (BELLE) [18] 



where Mi^2,3 are real and positive, and their numerical values are given in Eq. (|27|) - (|29|) . In our 
model, the Yukawa matrix is, 



/ 2.69 X io-6e-(0-695)i 5.92 x io-5e"(2-75)i 6.54 X io-4e-(i-68)i \ 



1.44 X l0-4e(i-54)» 1.73 X io-3e-(o i76)i g gi X io-3e-{i-32)i 
\^ 2.18 X io-3e(0-'^37)i 0.0213e(°-°°'^^)* 0.618 



(40) 



The non- vanishing off-diagonal matrix elements in ((5m~)jj induces lepton flavor violating processes 
mediated by the superpartners of the neutrinos through the one-loop diagram shown in Fig. ^ 

In Table ^ we summarize current status and future proposals of the experimental searches 
for lepton flavor violating decays. In the following subsections, we discuss each LFV process 
individually. 




Vi (5m^).j Vj 



J 



FIG. 1: The dominant diagram that contribute to the decay £i — > £jj at one loop, mediated by the neutralino 



XA and the sneutrinos v. The inserted mass term {Snij )ij is induced by the renormalization group evolution 



from the GUT scale to the RH neutrino mass scales. 
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67, r n'y and t —t 



The branching ratios for the decay of £j ^ £ j + 7 induced by the renormahzation group effects 
described above is given by [l2 1 

.3 _i 



^l^(3n.g + ^g)Ptan^/? 



fc=l,2,3 



Mr, ^' 



(41) 



Here a is the fine structure constant, Gp is the Fermi constant, and mg is the typical SUSY scalar 
mass which is given by, to a very good approximation [2^, 



1 



m 



/2) 



where M112 is the universal gaugino mass. In our model, |5(m|)j,| is given by. 



(42) 



\5{rr?j 



\l-{?,ml + Al)\ 

( * 3.41 X 10""^ 0.0098 ^ 

3.41 X 10"^ * 0.0962 

y 0.0098 0.0962 * j 



(43) 



for i ^ j- Thus the following relation is predicted, 



Br{iJ, 67) < Br{T — > 67) < Br(T H'y) . 



(44) 



Similar relation was observed in Ref. 



2ll | in which symmetric mass matrices with four texture zeros 
are utilized. We also note that the value for tan /? is 10, thus there is no tan /? enhancement in our 
predictions. 

Currently the most stringent experimental bound on the lepton flavor violating processes is on 
the decay fi — > 67. The prediction of our model for Br{fj, — > ej) is well below the most stringent 
bound up-to-date from MEGA at LANL 13]. In Fig. [2 the branching ratio of the decay /i — > 67 
as a function of the universal gaugino mass Mx/2 is shown for various scalar masses ^0 and rriQ. 
For large Aq and low tuq and M1/2, there is a large soft SUSY parameter space that give rise to 
predictions which can be probed by MEG at PSI and/or at J-PARC. In Fig. |31 the branching 
ratio of the decay r ^7 as a function of the universal gaugino mass M1/2 is shown for various 
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FIG. 2: The branching ratio of the decay ^ 67 as a function of the universal gaugino mass Afi/2 for 
various scalar masses Aq and mo. (SI): mo = Aq = 100 GeV; (Sll): mo = 100 GeV,Ao ^ 1 TeF; (S2): 
7^0 = Aq ^ 500 Ge^; (S3): mo = ^0 = 1 TeV. The dash line corresponds to the current experimental limit 
1.2 X 10^^^ from MEGA, while the solid line indicates the reach of a future experiment at J-PARC, 10~^^. 
The value of tan (3 in our model is tan j3 = 10. 
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FIG. 3: The branching ratio of the decay r — > /i7 as a function of the universal gaugino mass M1/2 for 
various scalar masses Aq and mo. (SI): toq = ^0 = 100 GeV; (Sll): mo = 100 GeV,Ao = 1 TeV; (S2): 
"m-o = j4o = 500 GeV; (S3): mo = Ao = 1 TeV. The dash line corresponds to the current experimental limit 
3.1 X 10^^ from BELLE, while the solid line indicates the reach of a future experiment at BELLE, 10~^. 
The value of tan/3 in om- model is tan/3 = 10. 
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scalar masses Aq and ruQ. For Aq ~ 0{1 TeV) and tuq and M1/2 both of order 0(100 GeV), the 
prediction of our model on r ^ fi'j may be tested at BELLE in the future. In Fig. [IJ the branching 
ratio of the decay r — > 67 as a function of the universal gaugino mass M1/2 is shown for various 
scalar mass Aq and niQ. For the SUSY parameter space we consider, the prediction for Br{T — > 67) 
is at least four orders of magnitudes below the current experimental upper bound. 
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FIG. 4: The branching ratio of the decay t ^ 67 as a function of the universal gaugino mass M1/2 for 
various scalar masses Aq and mo. (SI): niQ = Aq = 100 GeV; (Sll): mo = 100 GeV,Ao 1 TeV; (S2): 
TTT-o ^ Aq ~ 500 GeV; (S3): toq — Aq = I TeV. The dash line corresponds to the current upper bound, 
3.6 X 10^^, from BELLE. The value of tan/3 in our model is tan/3 = 10. 



We comment that, in models with lop-sided textures |22], the maximal mixing angle observed 
in the atmospheric neutrino sector is due to a large (23) mixing in the charged lepton sector. As a 
result, the off-diagonal elements in (23) sector of Oe^ are of order 0(1), which in turn gives rise to 
an enhancement in the decay branching ratios. In order to satisfy the current experimental upper 
bound, some new mechanism must be in place to suppress the decay rate of ^ C'j in models with 
lop-sided textures [2^. In our model which utilizes symmetric textures, as large leptonic mixing 
in our model is a result of the seesaw mechanism, all off-diagonal matrix elements in Y^, Oe^ and 
Or are much smaller than unity, leading to a much smaller branching ratio for /i — > 67 than that 
predicted in models with lop-sided textures. Yet our prediction is large enough to be probed by 
the next generation of experiments within a few years. 



B. /i 3e 



For the process /i 3e, as penguium diagrams are the dominant contributions, the branching 



'j ij'ij has similar structure as that of the decay 



ratio of the decay £• 

□ 

good approximation, the relation between these two processes reads |12l | 



' • 7. To a very 



Brii. 



s ^3 ^3 



a 
8^ 



3 ^2me/ 9 



where mi- is the i— th generation lepton mass. For the decay fj, — > 3e, we thus have 

Br{i^ ^ 3e) ~ 7 X 10"^Sr(/u ej) . 



(45) 



(46) 
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FIG. 5: The branching ratio of the decay fj,~ — > e~e"*"e~ as a function of the universal gaugino mass M1/2 
for various scalar masses Aq and tuq. (SI): ttiq = Aq ^ 100 GeV: (Sll): mo = 100 GeV, Aq = 1 TeV; (S2): 
niQ = ^0 = 500 GeV; (S3): mo — Aq ^ I TeV. The dash line corresponds to the current experimental limit 
1.0 X 10-12 from SINDRUM. The value of tan /3 in our model is tan /3 = 10. 

In Fig. [21 the branching ratio of the decay /i ^ 3e as a function of the universal gaugino mass 
Mx/2 is shown for various scalar mass Aq and tjiq. As the current experimental upper bound and 
the reach of the next phase of experiment at BELLE are still quite high, the prediction for /j, ^ 3e 
in our model can not be tested, even with a high value of the scalar mass, Aq = 1 TeV. 

C. /i-e Conversion 



Similar to the case of ;U — > 3e, the branching ratio for muon-electron conversion is also related to 
the branching ratio of the decay /x — > 67 as long as tan /3 is not too small. In the re gion tan/3 > 1, 
the relation between these two processes is given by, to a very good approximation 

Br{fi — > e) 



^^^^^.16aX,Z|F(,^)|^ 



(47) 



where -^e// is the effective charge of the nucleon, Z is the proton number and F{q^) is the nuclear 
form factor at momentum transfer q. For ^Ti, the conversion rate is 



Br[n e; \lTi) ~ 6 x 10"^Br(/x 67) , 



(48) 



where -^e// = 17-6 and F{q^ = —rnjj) = 0.54 have been used. In Fig. El the branching ratio of the 
decay — > e in ^Ti as a function of the universal gaugino mass Mi 12 is shown for various scalar 
.a. A. and F„ low ™K,es of and M.,„ the. . a v., l.ge sot. SUSY — e, spaoe 
that give rise to prediction for /u — e conversion rate that is sensitive to MECO |171 | at BNL and 
the proposal at J-PARC. 
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FIG. 6: The branching ratio of the decay fi~ — > e~ in IfTi- as a function of the universal gaugino mass M1/2 
for various scalar masses Aq and mo. (SI): tuq ^ Aq ^ 100 GeV; (Sll): niQ = 100 GeV, ^0 = 1 TeV; (S2): 
Too = Aq = 500 GeV; (S3): mo = Ao = 1 TeV. The dash line corresponds to the current experimental limit 
6.1 X 10^^'^ from SINDRUM II, while the solid line indicates the reach of a future experiment at J-PARC, 
2Q-i8^ The value of tan/3 in our model is tan/3 ~ 10. 

IV. BARYOGENESIS A LA SOFT LEPTOGENESIS 



It is well known that the CP violation in the quark sector is too small to explain the observed 
baryon asymmetry of the Universe (BAU), expressed in terms of the ratio of the baryon number 
to entropy 



rib 



(0.87 ±0.04) X 10 



-10 



(49) 



derived from CMB and nucleosynthesis measurements. In leptogenesis, leptonic CP violating 
phases are used to produce asymmetry in leptonic number which then is converted into baryon 
asymmetry by the electroweak non-perturbative effects due to sphalerons. There are two ways of 
producing lepton number asymmetry: (i) Standard leptogenesis (STDL) j^l and (ii) Soft leptoge- 
nesis (SFTL) 22, 3- 

In STDL scenario, the primordial leptonic asymmetry is generated by the decay of the heavy 
right-handed Majorana neutrinos and their scalar partners, mediated by the Yukawa interactions 
in the superpotential. In our model, the large hierarchy among the three heavy neutrinos leads 
to a very small CP asymmetry, which is of the order of 0(10^^). In addition, the low value for 
the mass of the lightest RH neutrino. Mi = 1.09 x 10^ GeV, leads to an extremely large wash-out 
effect. Due to these reasons, the prediction in our model for the baryonic asymmetry utilizing the 
standard leptogenesis is of the order of 0(10^^^), which is four orders of magnitude below the value 
derived from experimental observations. 
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SFTL utilizes the soft SUSY breaking sector, and the asymmetry in the lepton number is 
generated in the decay of the superpartner of the RH neutrinos 



26 



23], as opposed to the hghtest 



RH neutrino in the case of STDL. UnUke in STDL where the Yukawa sector is responsible for 
the required CP violation and lepton number violation, in the scenario of SFTL, the CP violation 
and lepton number violation trace their origins to SUSY breaking. As a result, it allows a much 
lower bound on the mass of the lightest RH neutrino, Mi, compared to that in STDL. In fact, it 
has been shown very recently that in contrast to the STDL scenario in which Mi > lO^GeV is 
typically required to have sufficient baryonic asymmetry i29|, SFTL can only work in the region 

. As a result, the problem of the gravitino over-production 31] may be 

avoided. 

For SFTL, the relevant soft SUSY Lagrangian that involves lightest RH sneutrinos i^r-^ is the 
following, 

- Csoft = i^BMiUR^UR^ + AyiiLiUR^Hu + h.c.) 

+fr?V'^j^VR^ . (50) 

This soft SUSY Lagrangian and the superpotential that involves the lightest RH neutrino, A^i, 

W = MiNiNi + yuUNiHu (51) 

give rise to the following interactions 

-C^ = VR, (MiYUlK + yx^JL + AyxMu) 

+h.c. , (52) 

and mass terms (to leading order in soft SUSY breaking terms), 

-C^ = {mIv\Vr^ + ^BMiVrJ)r^) + h.c. . (53) 

Diagonalization of the mass matrix M with the two states vr-^ and z?|j^ leads to eigenstates 
and with masses, 

M±=.Mi(l±i^), (54) 

where the leading order term Mi is the F-term contribution from the superpotential (RH neutrino 
mass term) and the mass difference between the two mass eigenstates and is induced by the 
SUSY breaking B term. The time evolution of the T^Ri-i^'^r^ system is governed by the Schrodinger 
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equation, 

where the Hamiltonian Tl is given by 26, 2^ . 

n = M-^£/ (56) 

M = ( Ml, (57) 

(58) 

For the decay of the hghtest RH sneutrino, the total decay width Fi is given by, in the basis 
defined in Eq. p7() where both the charged lepton mass matrix and the RH neutrino mass matrix 
are diagonal, 

Tl = -^{yuylhiMi = 0.374 GeV . (59) 
47r 

The eigenstates of the Hamiltonian 7i are -/Vj_ = pN it qN^ , where |pp + |gp = 1. The ratio q/p is 
given in terms of M and F as. 




^ 2M 



12 '"^12 



p ) 2Ml2 — i£^12 

-i+m^). (60) 

in the limit £/i2 <^ Aii2- Similar to the — system, the source of CP violation in the lepton 
number asymmetry considered here is due to the CP violation in the mixing which occurs when 
the two neutral mass eigenstates (A^+, -/V-), are different from the interaction eigenstates, (A^+, 
N'_)- Therefore CP violation in mixing is present as long as the quantity \q/p\ 7^ 1, which requires 

For this to occur, SUSY breaking, i.e. non-vanishing A and B, is required. As the relative 
phase between the parameters A and B can be rotated away by an [/(l)/j-rotaion, without loss of 
generality we assume from now on that the physical phase that remains is solely coming from the 
tri-linear coupling A. 

The total lepton number asymmetry integrated over time, e, is defined as the ratio of difference 
to the sum of the decay widths F for ur-^ and t'jj^ into final states of the slepton doublet L and the 
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Higgs doublet or the lepton doublet L and the higgsino H or their conjugates, 



(62) 



where final states f = {L H), {L H) have lepton number +1, and / denotes their conjugate, 
H^), {L H), which have lepton number —1. After carrying out the time integration, the total 



CP asymmetry is 



26 



3, 

ATiB \Im{A) 



where the additional factor 5b-f takes into account the thermal effects due to the difference 



between the occupation numbers of bosons and fermions 
asymmetry is 



26 



jaiic 

3, 



3- 



The final result for the baryon 



riB , 
— ~ -cd^ e K 
s 



-1.48 X lO^^e K 



[lA8xlO-^)(^^^ R5B-FK (64) 



where in the first line is the density of the lightest sneutrino in equilibrium in units of entropy 
density, and is given by, dp^ = 45C(3)/(7r'^5*); the factor c = (SNp + ANh)/{22Nf + 13iV//) 
characterizes the amount oi B — L asymmetry being converted into the baryon asymmetry Yb, 
with Np and Nh being the number of families and the SU{2) Higgs doublets, respectively. For 
the MSSM particle spectrum, {Nf-,Nh) = (3,2). The parameter k is the dilution factor which 
characterizes the wash-out effects due to the inverse decays and lepton number violating scattering 
processes together with the time evolution of the system. It is obtained by solving the Boltzmann 
equations for the system. An approximation is given by [3 

10^ < r : K = (0.1r)i/2e-(|)(0-H^/* (65) 
10<r<10^: K = 0.3/(r(lnr)°-6) (66) 



< r < 10 : K= l/(2Vr2 +9) . (67) 

where r is defined as 



(1.7) (32^) Ml 

with Mpi being the Planck scale taken to be 1.2 x 10^^ GeV . We have r = 183 and correspondingly 
K, = 0.00061 in our model. The parameter R is defined as the ratio, 
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which gives a value equal to one when the resonance condition, Fi = 2\B\, is satisfied, leading 
to maximal CP asymmetry. As Ti is of the order of 0(0.1 — 1) GeV, to satisfy the resonance 
condition, a small value for B <^ fh \s thus needed. Such a small value of B can be generated 
by some dynamical relaxation mechanisms in which B vanishes in the leading order. A small 
value of i? ~ fh? /Ml is then generated by an operator J d'^OZZ'^N'^/M^^ in the Kahler potential, 
where Z is the SUSY breaking spurion field, Z = 9'^ fhMpi [2^. In our model, with the parameter 
B' = \/BMi having the size of the natural SUSY breaking scale Vm? ~ 0(1 TeV), a small value 
for B required by the resonance condition S ~ Fi ~ 0(0.1 GeV) can thus be obtained. 

Fig. Q shows the ratio ii as a function of B'. For the specific value of the decay width Fi 
predicted in our model, the resonance occurs at around B' ~ 1.4 TeV. In Fig. |HJ the region on 
the Im{A) versus B' plane that gives rise to an amount of baryon asymmetry consistent with the 
value derived from observation, ub/s = (0.87 it 0.04) x 10^^^, is shown. The required value for 
B' near the resonance is around 800 GeV — 2 TeV, and the required value for |/m(A)| is around 
(1 — 2) TeV. At the resonance B', the value for |/m(74)| can be as low as 1 TeV to generate 
sufficient amount of baryon asymmetry. In Fig. IHl we show the predictions for the asymmetry, 
ub/s, as a function of B' for different values of Im{A). In the numerical analyses presented in 
Fig. ISl and IHl we assume 6b-f = 1- We note that even if an additional suppresion 6b-f ~ 0.1 
is present, with a value of Im(A) ~ 10 TeV at the resonance our model can still account for the 
observed BAU. 
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FIG. 7: The ratio i? as a function of B' . The resonance occurs at around B' ^ l.A TeV . 
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FIG. 8: The parameter space on the Im{A) versus B' plane that gives rise to an amount of baryon 
asymmetry consistent with the value derived from observations, ub/s — (0.87 ± 0.04) x 10^^*^, is the 
region bounded by these two curves. The upper curve corresponds to the upper bound from observation, 
ub/s = 0.91 X 10^^°, while the lower curve corresponds to the lower bound, ns/s = 0.83 x 10~^°. 
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FIG. 9: The prediction for ub/s as a function of B' for |/to(^)| = 10 TeV (circles), 5 TeV (triangles) and 
1 TeV (squares). 

V. CONCLUSION 

We have shown in this paper that, in contrast to the predictions of models with lop-sided 
textures, the predictions for LFV decays are well below the current experimental bounds. This is 
demonstrated in a model based on SUSY SO(IO) with symmetric mass textures which give rise to 
predictions for all fermion masses and mixing angles, including those in the neutrino sector, that 
are in good agreement with experimental data within la. The predictions of our model for LFV 
processes, ii — > Ij^, /x — e conversion as well as /i — > 3e, are well below the most stringent bounds 
up-to date. Our predictions for many processes are within the reach of the next generation of 




lm(A) = 10 VeV cT 
lm(A) = 5TeV a. 
<^ V lm(A) = 1 TeV ❖ 

f> '>c 
c \ 

«> X„ 

- « A-" "'^-A 

« A- "^A^A,, 



C A" "A^A. 
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LFV searches. This is especiahy true for /x — e conversion and fj, 67. We have also investigated 
the possibiUty of baryogenesis resulting from soft leptogenesis. Our model predicts Mi < lO^GeV 
which is the required condition for this mechanism to work. With the soft SUSY masses assuming 
their natural values, B' ~ 1.4 TeV and Im{A) ~ 1 TeV, we find that our model can indeed 
accommodate the observed baryon asymmetry of the Universe. 

Acknowledgments 

M-CC and KTM arc supported, in part, by the U.S. Department of Energy under Grant No. DE- 
AC02-98CH10886 and DE-FG03-95ER40892, respectively. M-CC would also like to acknowledge 
Aspen Center for Physics, where part of this work was done, for its hospitality and for providing 
a very stimulating atmosphere. 



[1] M.-C. Chen and K. T. Mahanthappa, Int. J. Mod. Phys. A18, 5819 (2003). 

[2] M.-C. Chen and K. T. Mahanthappa, Phys. Rev. D62, 113007 (2000); ibtd. D65, 053010 (2002); ibtd. 
D68, 017301 (2003). 

[3] R. Barbieri, L. J. Hall, S. Raby, and A. Romanino, Nucl. Phys. B493, 3 (1997). 
[4] S. Eidelnian et al. (Particle Data Group), Phys. Lett. B592, 1 (2004). 

[5] H. Fusaoka and Y. Koide, Phys. Rev. D57, 3986 (1998); A. Hocker, H. Lacker, S. Laplace, and 

F. Le Diberder, Eur. Phys. J. C21, 225 (2001). 
[6] J. Charles ct al. (CKMfitter Group), hep-ph/0406184; M. Bona ct al. (UTfit), hep-ph/0408079. 
[7] L. J. HaU, R. Rattazzi, and U. Sarid, Phys. Rev. D50, 7048 (1994). 

[8] M. Mahoni, T. Schwetz, M. A. Tortola, and J. W. F. Valle, New J. Phys. 6, 122 (2004), hep-ph/0405172. 
[9] J. N. BahcaU, M. C. Gonzalez-Garcia, and C. Pena-Garay, hep-ph/0406294. 
[10] M. V. Diwan et al., Phys. Rev. D68, 012002 (2003). 

[11] For a recent review, see A. Masiero, S. K. Vempati, and O. Vives, hep-ph/0407325. 

[12] J. Hisano, T. Moroi, K. Tobe, and M. Yamaguchi, Phys. Rev. D53, 2442 (1996); J. Hisano and 

D. Nomura, Phys. Rev. D59, 116005 (1999). 
[13] M. L. Brooks et al. (MEGA), Phys. Rev. Lett. 83, 1521 (1999). 

[14] L. M. Barkov et al., Research Proposal to PSI, http:/ /www.icepp.s.u-tokyo.ac.jp/meg (1999). 
[15] U. Bellgardt et al. (SINDRUM), Nucl. Phys. B299, 1 (1988). 

[16] P. Wints (SINDRUM II), in Proceedings of the First International Symposium on Lepton and Baryon 
Number Violation, ed. H. V. Klapdor-Kleingrothaus and I.V. Krivosheina, p. 534 (1998), Institute of 
Physics Publishing, Bristol and Philadelphia. 



21 



[17] M. Bachman et al. (MECO), Proposal to BNL, http://meco.ps.uci.edu (1997). 

[18] K. Inami (BELLE), Talk presented at the 19th International Workshop on Weak Interactions and 

Neutrinos (WIN'03), Lake Geneva, Wisconsin, (2003). 
[19] K. Abe et al. (Belle), Phys. Rev. Lett. 92, 171802 (2004). 

[20] S. T. Petcov, S. Profumo, Y. Takanishi, and C. E. Yaguna, Nucl. Phys. B676, 453 (2004). 

[21] M. Bando, S. Kaneko, M. Obara, and M. Tanimoto, hep-ph/0405071. 

[22] C. H. Albright, K. S. Babu, and S. M. Barr, Phys. Rev. Lett. 81, 1167 (1998). 

[23] S. M. Barr, Phys. Lett. B578, 394 (2004); T. Blazek and S. F. King, Nucl. Phys. B662, 359 (2003). 

[24] C. L. Bennett et al., Astrophys. J. Suppl. 148, 1 (2003). 

[25] M. Fukugita and T. Yanagida, Phys. Lett. B174, 45 (1986); M. A. Luty, Phys. Rev. D45, 455 (1992); 

W. BuchmuUer and M. Plumacher, Phys. Lett. B389, 73 (1996). 
[26] Y. Grossman, T. Kashti, Y. Nir, and E. Roulet, Phys. Rev. Lett. 91, 251801 (2003). 
[27] G. D'Ambrosio, G. F. Giudice, and M. Raidal, Phys. Lett. B575, 75 (2003). 

[28] L. Boubekeur, hep-ph/0208003; L. Boubekeur, T. Hambye and G. Senjanovic, Phys. Rev. Lett. 93, 

111601 (2004), hep-ph/0404038. 
[29] W. BuchmuUer, P. Di Bari, and M. Plumacher, Phys. Lett. B547, 128 (2002); G. F. Giudice, A. Notari, 

M. Raidal, A. Riotto, and A. Strumia, Nucl. Phys. B685, 89 (2004). 
[30] Y. Grossman, T. Kashti, Y. Nir, and E. Roulet (2004), hep-ph/0407063. 

[31] M. Y. Khlopov and A. D. Linde, Phys. Lett. B 138, 265 (1984); B. A. Campbell, S. Davidson, and K. A. 

Olive, Nucl. Phys. B399, 111 (1993); S. Sarkar, hep-ph/9510369; G. G. Ross and S. Sarkar, Nucl. Phys. 

B461, 597 (1996); M. Y. Khlopov, Cosmopartide physics, 577p, World Scientific, Singapore (1999). 
[32] L. Covi, N. Rius, E. Roulet, and F. Vissani, Phys. Rev. D57, 93 (1998). 

[33] E. W. Kolb and M. S. Turner, The Early Universe, (1990), Redwood City, USA: Addison- Wesley 547 

p. (Frontiers in physics, 69); H. B. Nielsen and Y. Takanishi, Phys. Lett. B507, 241 (2001). 
[34] M. Yamaguchi and K. Yoshioka, Phys. Lett. B543, 189 (2002). 



22 



